Can modified gravity explain accelerated cosmic expansion? 
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l/"") \ Abstract 

, We show that the recently suggested explanations of cosmic acceleration by the mod- 

^ ' ification of gravity at small curvature suffer violent instabilities and strongly disagree 

, with the known properties of gravitational interactions. 
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It seems established that at the present epoch the universe expands with acceleration. 
£T) \ It follows directly from the observation of high red-shift supernovae [1] and indirectly from 

the measurements of angular fluctuations of cosmic microwave background fluctuations 
(CMBR) [2]. The latter shows that the total mass/energy density of the universe is very 
close to the critical one (f2 = 1), while the observations of the universe large scale structure 
indicate that normal gravitating (but invisible) matter can contribute only 30% into the 
total one. Thus one concludes that the remaining 70% is some mysterious agent that creates 
the cosmological acceleration. The simplest suggestion is that the source of this acceleration 
is the vacuum energy (cosmological constant). However such an explanation immediately 
meets two serious problems. First, vacuum energy remains constant in the course of cosmic 
^ ■ expansion (modulo possible phase transitions which could discontinuously change the value 

°f Pvac by a very large amount). On the other hand, "normal" cosmological energy density 
scales with time as p m ~ m 2 Pl /t 2 and it is quite strange coincidence that p vac and p m are of 
the same order of magnitude just today. Second and much more serious, there are known 
contributions into vacuum energy which are 50-100 orders of magnitude larger than the 
cosmologically allowed value. At the present day there is no workable scenario to solve 
the second problem, while the first one may be solved if instead of a constant p vac a new 
massless or very light field is introduced which evolves in the course of the expansion in such 
a way that its energy density today more or less naturally approaches the energy density of 
normal matter [3]. 

Still such field remains mysterious. It would be in much better shape if it simultaneously 
solved the problem of almost complete elimination of gigantic contributions into vacuum 
energy by some adjustment mechanism but, alas, it does not happen or, better to say, it 
is unknown how it can be realized. To avoid an introduction of additional fields (plurality 
should not be posited without necessity" , W. Ockham) it was suggested recently that gravity 
itself, if properly modified, could create cosmological acceleration [4]. To this end the model 
with the following action was considered: 

A= /<fV=9 [/(«) + £«] , (1) 
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where R is the scalar curvature and Cm is the matter Lagrangian. In the usual Einstein 
gravity the function f(R) has the form f(R) = Rm 2 pjl&it where mpi is the Planck mass 
(~ 10 19 GeV). In ref. [4] a power law modification f(R) ~ R n with n = —1,3/2 was 
considered. 

A detailed examination of the n = —1 case was performed in the recent paper [5]. The 
action was taken in the form: 

A = it J d ' x ^ ( R ~ 7f) " / * x ^ Cm ' (2) 

and the corresponding equation of motion reads: 

(4 \ -j / 4 \ /IN /IN 8ir 

1 + ^ - ~ 2 (l - % ) «M + (si) - «*V(.V» (^) = • (3) 

where is the matter energy- momentum tensor. As argued in ref. [5], with the parameter 
H chosen of the order of the inverse universe age, fj,^ 1 ~ 10 18 sec ~ (10 _33 eV) _1 , this 
equation, applied to cosmology, describes the universe acceleration in the present epoch, 
while the additional terms were not essential at earlier times. 

However an addition of non-linear in curvature terms into the action integral is not 
an innocent step. It gives rise to higher order equations of motion which possess quite 
unpleasant pathological behavior. They may break unitarity and lead to ghosts or tachyons. 
And if breaking of the usual gravity at high curvatures may be tolerated - who knows what 
happens there - but breaking of gravity at low curvature in weak field regime immediately 
leads to contradiction with well established facts. To demonstrate that let us consider 
equation of motion for curvature scalar. It can be obtained from eq. (3) by taking trace 
over a and /?: 

2 (D a R) (D a R) R 4 R 2 TR* 
DR ~ 3 R + 6V"^" = "^ r ' () 

where D is the covariant derivative and T = 8ttT" /m Pl > 0. 

Let us apply this equation for the gravitational field of a normal celestial body, as e.g. 
the Earth or the Sun or any other smaller gravitating object with weak gravity. In this case 
the metric can be approximately taken as a flat one, so D 2 = d 2 — A and (D a R)(D a R) = 
(dfR) 2 — (djR) 2 . We will look for a perturbative solution inside some spatially finite matter 
distribution. In the lowest order the curvature is algebraically expressed through the trace 
of the energy-momentum tensor of matter, Ro = —T. This is the standard result of General 
Relativity (GR). Outside matter the solution is Rq = but now this is not exact but true 
only in the lowest order approximation. One can check both numerically and analytically 
that stationary solution outside matter sphere is very quickly vanishing if fi 4 > 0. So at least 
in this case the curvature scalar behaves similarly to the usual one. This is in agreement with 
the analysis of stationary solutions of ref. [6] where it is argued that the modified gravity 
theories agree with the Newtonian limit of the standard gravity for sufficiently small [i. 1 

However for negative fi 4 the solution is rising outside the matter sphere and strongly 
deviates from the standard General Relativity. Even for /j, 4 > 0, the curvature scalar, 
R, outside the source would be non-vanishing and proportional to the curvature on the 

1 In Ref. [6], however, it was found that the condition for the correct Newtonian limit is very marginally 
satisfied for the simple 1/R gravity. Thus, it may be likely that the deviation from the Newtonian gravity 
is observed even in vacuum solutions. 
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boundary of the gravitating body but, as we mentioned above, it tends to zero at very 
small distances from the boundary. 

Let us consider now the internal solution with time dependent matter density. The first 
order correction to the curvature, R = —T + Ri satisfies the equation: 



Ri - AR 1 -^R 1 + ^d j R 1 + R 1 

T 2 3D a TD a T 
= D 2 T + — , 



rr i c> T 2 — (djT) 2 T 3 
1+6 T 2 6/x 4 



(5) 



where the deviation of the metric from the Minkowski one may be accounted for in the 
r.h.s. by the covariant derivatives (but this is not essential for our conclusion). 
The last term in the l.h.s. has a huge coefficient: 

T 3 /6/ ~ (l(r 26 fc 2 sec)- 2 ( -^-V) 3 , (6) 

\g cm.- 6 J 

where p m is the mass density of the gravitating body, the constant k is assumed to be of the 
order of unity, fi = k/t u and t u ~ 3 ■ 10 17 sec is the cosmological time scale. This coefficient 
is much larger than T since 

T ~ (10 3 sec)- 2 ( • (7) 
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Thus, the T 3 /6/i 4 dominates the coefficient in front of R\ term in Eq. (6) and leads to 
strong instability. (Notice that —T 3 /6fi 4 < 0.) The characteristic time of the instability 
rise is about 1CT 26 sec. Perturbation with the wave length larger than one tenth of the 
proton Compton wave length should be unstable. This result shows that the curvature 
inside matter sphere should quickly rise and reach very high values (strong gravity) till 
some non-perturbative effects would terminate this rise. This is the same kind of instability 
that leads to the universe acceleration but now it is very unfavorable for the model and 
forbids such a modification of gravity. 

It was argued in a recent work [7] that the considered above modification of gravity 
follows from certain compactification schemes of string/M-theory. If so, such theories must 
be excluded because of the instability found in the present paper. 

The analysis performed above is valid for a particular case of 1/R terms in the action 
integral. A more general Lagrangian with arbitrary powers of R, both negative and positive, 
is presented in the recent work [8]. We have not made complete analysis of all possible cases 
discussed in the literature but it seems suggestive that all modifications of gravity at small 
curvature described by higher order differential equations may suffer from serious instability 
problems. 

There are other earlier works on modifications of gravity [9]- [13] which might lead to 
accelerated cosmological expansion. However it is impossible to apply the similar analysis 
to some of them because the theory is not sufficiently specified - only a modification of 
the Friedmann equations is suggested but no modified action which would allow to derive 
equations in general case was presented [9]-[ll]. Another group of works is based on modifi- 
cation of gravity due to higher dimensions [12]-[13]. Probably these models do not manifest 
instability discussed here, though only cosmological solutions have been considered in the 
quoted papers. In ref. [14] it is argued that modification of gravity at large scales to in- 
duce cosmic acceleration can lead to a modification of gravitational interactions at smaller 
scales but since the authors did not specify their model it is hard to make any definitive 
conclusions. 
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To summarize, we have shown that an attempt to explain the observed cosmic accel- 
eration by the concrete modification of gravitational action adding to the latter the term 
~ 1/R is unacceptable because such term would lead to a strong temporal instability result- 
ing in a dramatic change of gravitational fields of any gravitating bodies. We cannot make 
any conclusion about some other versions of the modified gravity theory which could lead 
to cosmic acceleration because the latter are not sufficiently well defined to permit such an 
analysis. 
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